IPPP/06/64 
DCPT/06/128 
October 2006 



Improved analysis of moments of in neutrino-nucleon 
scattering using the Bernstein polynomial method 

Paul M. Brooks* and C.J. Maxwell^ 



Institute for Particle Physics Phenomenology, University of Durham, 
South Road, Durham, DHl 3LE, UK. 



Abstract 

We use recently calculated next-to-next-to-leading order (NNLO) anomalous dimen- 
sion coefBcients for the moments of the xFy, structure function in vN scattering, together 
with the corresponding three-loop Wilson coefficients, to obtain improved QCD predic- 
tions for both odd and even moments of F^. To investigate the issue of rcnormalization 
scheme dependence, the Complete Rcnormalization Group Improvement (CORGI) ap- 
proach is used, in which all dependence on rcnormalization and factorization scales is 
avoided by a complete resummation of RG-prcdictable scale logarithms. We also con- 
sider predictions using the method of effective charges, and compare with the standard 
'physical scale' choice. The Bernstein polynomial method is used to construct experimen- 
tal moments (from the xF^ data of the CCFR collaboration) that are insensitive to the 
value of xF^ in the region of x which is inaccessible experimentally. Direct fits for 
{as{Mzy) are then performed. The CORGI fits including target mass corrections give a 
value as{Mz) = O.IISOIqqqJq, consistent with the world average. The effective charge 
and physical scale fits give slightly smaller values, which are still consistent within the 
errors. 
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1 Introduction 



The measurements of the CCFR cohaboration provide a precise determination of the non- 
singlet deep inelastic scattering (DIS) structure functions of neutrinos and anti-neutrinos on 
nucleons, xF3{x, Q^) [1]. Recently the NuTev collaboration have also published measurements 
of F2{x, Q"^) and F^{x^ Q^) [2]. The Q^-dependence of moments of structure functions can be 
predicted in perturbative QCD, and fits to the data can be used to infer Ams (or equivalently 
as{Mz))- In doing this, the principal difficulty is that there are upper and lower limits on 
the experimentally accessible range of x at low and high Q^, respectively. The moments are 
potentially sensitive to this missing information, and this propagates into an additional level 
of uncertainty in the resultant prediction of Ams. 

An approach which has been applied in the past is to use Bernstein polynomials, 
which are peaked in a rather limited x-range, to construct linear combinations of moments 
which are insensitive to the missing x regions [3l Hj. The analysis of Ref. [3] chose, as is 
customary, to work in the MS scheme and set both renormalization and factorization scales 
to Q (physical scale (PS) choice). This was extended in Ref. [5] to consider predictions for F3 
obtained in the 'complete renormalization group improved' (CORGI) approach [6] in which 
all dependence on the renormalization scale /u and the factorization scale M is eliminated by 
an all-orders resummation of RG-predictable scale logarithms. 

The analyses of Refs. [31 14] and [5] used the then state-of-the-art three-loop (NNLO) 
results for the anomalous dimension and coefficient function, which were restricted to a subset 
of odd moments n = 1, 3, 5, . . . , 13 [7]. Recent progress has yielded NNLO results for these 
quantities for any value of n [HI [9] . Consequently the set of Bernstein moments used in the 
fits can now be greatly extended. The Bernstein polynomials defined in Refs. [3l Hj were 
linear combinations of odd moments, but the new results of Refs. [8l[9] mean that, by a slight 
redefinition of the polynomials, even moments can now also be studied. 

In this paper we intend to perform such an extended analysis. We shall fit the CCFR 
data [1] to PS and CORGI NNLO QCD predictions, and will in addition compare with the 
predictions in the closely-related method of effective charges approach (EC) [10]. Target mass 
corrections and higher twist effects will also be considered. We shall also compare our results 
to those obtained using a fitting technique based on Jacobi Polynomials [111 [T2| [T3] . 

The plan of the paper is to give a brief review of the factorization and renormal- 
ization scheme dependence of structure function moments in Section 2. In Section 3 we 
discuss the CORGI and effective charge approaches for leptoproduction moments. We take 
this opportunity to correct an error in the expression for the NNLO CORGI result for the 
scheme invariant X2 derived in Ref. [5] . Section 4 will contain a description of the Bernstein 
polynomial averages to be employed in the fits. We shall show how to modify the definition 
of the polynomials to accommodate both odd and even moments. We then constrain the 
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set of acceptable Bernstein moments to be used in the fits by comparing how four different 
methods of extrapolation (to obtain xF^ on the full x-range) differ; this enables us to define 
a 'modelling error' to be combined with the other sources of error in our analysis. In section 
5 we give details of the fitting procedure and in section 6 we present the results of the fits 
to the PS, CORGI and EC predictions for the moments, and consider how the fits change 
if target mass corrections and higher twist corrections are included. Section 7 contains a 
discussion and conclusions. 

2 Factorization and renormalization scheme dependence of 
the moments 

The moments we are concerned with in this paper are those derived from in (anti)neutrino- 
nucleon scattering. They are defined as: 

A^3-^(n;Q2) = [\x x^"'F^^ ix,Q^). (1) 
Jo 

These moments can be factorized in the following form, 

Mfin;Q^) = {N\On, ^s{M)\N) C^^Hq, M, ,2, a{fi)), (2) 
where {N\On.NsiM)\N) is the non-singlet (NS) operator matrix element of nucleon states and 

(3) 

Cn {Q, M, ^,a{^)) is the coefficient function. Here a = a^/vr is the RG-improved coupling. 
The operator matrix element is factorized at the scale M into a non-perturbative component 
and a perturbative expression, written in terms of the coupling evaluated at the factorization 
scale, a = a{M). The factorization scale dependence is governed by the anomalous dimension 
equation, 

M^{N\On,ns{M)\N) = {N\On,^s{M)\N)^n,^s{a). (3) 

Here Ns(fl) is the anomalous dimension of the moment. It has the following perturbative 
expansion, 

7n, Ns(a) = -d{n)a - di{n)a'^ - d2{n)a^ - d^{n)a^ - . . . , (4) 

where d{n) is factorization scheme invariant, and the higher coefficients serve to label the 
factorization scheme dependence. The M-dependence of the coupling is governed by the 
beta-function equation, 

^^TTTT = /3(a) = -ha^{l + ca + C2a^ + 030?- + ...). (5) 
oM 

Here b = {2>'i — 2Nf )/Q and c = (153 — 19iVj;)/126 are renormalization scheme (RS) invariant. 
The higher coefficients serve to label the RS dependence. Together, these two equations 
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determine the perturbative behaviour of the operator matrix element. For the remainder of 
this paper we simphfy our notation by dropping the sub- and superscripts ^vN\ 'n', '(3)' and 
'NS', from the quantities in Eqs. ([2|) and ([3j). Also, although the coefficients di{n) in Eq. ([4]) 
are n-dependent, we also suppress this. 

A solution to Eq. ([3]) can be obtained in the form, 

{0(M)) = ^„exp|^ 1^..-^ 2^..|. (6) 

where 7*^*) and /J*^*^ denote the anomalous dimension and beta-function equations truncated 
after i terms. There is a distinct parallel between the above equation and the solution to 
the beta function equation. The second integral in Eq. Q is an infinite constant. We 
are free to choose any form we wish for this term, subject to the constraint that it must 
have the same singularity structure as the first integral. However, a particular choice for 
this constant corresponds to a particular definition of An- Consequently, An can be likened 
to the dimensional transmutation parameter. A, in that it defines the missing boundary 
condition in Eq. ([3]). An is actually a (set of) non-perturbative constant(s), generated by the 
factorization process, and they are factorization and renormalization scheme (FRS) invariant. 
Their precise values cannot be calculated within perturbation theory, and hence must be 
obtained by comparison with experimental data. 

The coefficient function C{Q, M, fj,,a{fi)), depends on both the renormalization and 
factorization scheme adopted, and it takes the form of an expansion in powers of the coupling 
evaluated at the renormalization scale, 

C{Q,M,fi,d{fi)) = l + rid + r2a^ + r3d^ + (7) 

where d = a{M = fj,). Using the above equation together with Eqs. ([2]) and ([6]), the moments 
can be written as [U [H] , 

/ \ d/b 

M{n;Q^) = A„ ( — ^ ) exp(T(a))(l + na + rgoV rga^ + • ■ •), (8) 
\l + caj 

where, 

r r di + {dic + d2- dc2)x + ((ig + cdg - Cgfi)^^ + • • •) 
l[a) = j dx 15 ^ . (9) 

Jo 0(1 + CX)(1 + CX + C2X^ + CsX'^ H ) 

The explicit M dependence of the coupling can be obtained by solving the following tran- 
scendental equation [H], 

1 1 



1 , ca M ^ 

— hem = bm— 

a 1 + ca A 







/3(x) /?(2)(x) 



(10) 



Equation ([8]) serves as a prototypical expression for the moments, from which CORGI, EC 
and PS predictions can be derived. 
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The self-consistency of perturbation theory means that the perturbative coefficient 
ri has a dependence on M and di. The higher coefficients also have a dependence on the 
parameters specifying the FRS, M, C2, C3, . . ., c^; di,d2, ■ ■ ., c^fe)- The explicit form of this 
FRS-dependence can be determined by demanding that on calculating the moments up to 
0(0^^) the partial derivative with respect to each FRS parameter is 0(0^"*"^) [151 IS] ■ The 
complete set of partial derivatives required to derive the FRS-dependence of ri, r2 and r-^ is, 
for the /x-dependence 

= 0, f^^=nb, f^— = 2r2b + bnc. (11) 

For the M-dependence we have 

^r^ dro 

M = d2 + dm + dr2 - dnL - 2diL + dL"^ - dcL, (12) 
dM 

where we have defined L = bln{M/fi). For the C2-dependence we have 



For the cs-dependence 

dri 
dc3 

For the di-dependence 



dri 
dc2 


dr2 

dC2 


d 

~ "26' 




drs 
dc2 


rid Ld 
^ ~^^~b' 


cd 2di 
^ 36 ~ "3^ 


- n- 


= 0, 


dro 

— - = 0) 


drs _ 
dc3 


d 
~Qb' 


dri 

Wi " 


1 dr2 
~V ddl 


c ri 
^ 2b~J 


L 


dr^i 
ddi 


cri (? r2 
'2b ~3b~T 


C2 

36 6 





(13) 



(14) 



(15) 



For the d2-dependence we have 



dri dr2 1 dr^ c 1 

5d^ = °' dd-2=-Yb' 9d^ = 36 + 26^'^- ^''^ 
Finally, for the da-dependence we have 

^ = ^ = ^ = -1 (17) 
dds ' 5d3 ' 9^3 36' ^ ' 

These results may now be integrated to obtain ri, r2 and r^. For ri we obtain 

ri = ^TM-^-Xo{Q), (18) 
6 
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where tm = hln ^M/Aj. Xq{Q) is an FRS invariant quantity, generated as a constant of 
integration. One can define an FRS invariant, non-universal scale parameter, A_a/), via the 
FRS invariant Xq{Q). Thus, 

= X„(Q)..l„(^). (19) 

For r2 we obtain 

/ 1 b \ 2 b ~ di dc2 cdi d2 „ 

r2 = o-TT3ki + 3^i^i + -7^i-^ + TTr: + ^-TTr+^2, (20) 




where we have defined. 



d d 2b 2bd 2b 2b 



fi = ri{M = fi) 



= \r,-^-Xo{Q). (21) 

Here X2 is another FRS-invariant constant of integration. Crucially X2 and higher invariants 
are independent of Q. Hence, the complete Q-dependence of the observable is generated by 
Xq{Q). Similarly, for r3 we obtain 

C2dc c^d cd2 ds 2df cdf c^di did2 02^1 
~ ~W ~ ~6b ^ ^ ~ 3b~ 3bd? ~ 2bd~ "sT ^ Ixf ^ ~W 
2d\ri dirf b'^rf rf bcrf d2ri 

2bdirifi b'^rirf brffi bcrifi 

d2 d2 d~^~d~ 

2br^r2 2d-ir2 , , 

+ -^ + ^^+rir2 + X^. (22) 

Again X3 is a Q-independent FRS-invariant constant of integration. Using Eq. ()20p X3 can 
be written in terms of ri, fi and the other FRS parameters. This also holds for the higher 
invariants. The results of Eqs. ([IT]) - (fT7|) and of Eqs. ([20]) and (p2]) replace, respectively, 
Eqs. (15) and (18) of Ref. [5] which contain several errors. The invariant X2 can be obtained 
from NNLO results for the anomalous dimension coefficients and coefficient function in any 
FRS. For instance if we make the customary choice of MS with Af = ^ = Q then ri = f\ and 
we obtain 




^ ^ di dc2 d\ cdi d2 



_• (23) 

MS 

In summary, through Eqs. (jlSp . (j20p and (|22p . we have determined the explicit FRS 
dependence of the coefficients r^. In doing so, we have generated a set of FRS invariant 
quantities Xi, the importance of which will become clear when we come to consider the 
CORGI form of the moments in the following section. 
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3 PS, CORGI and EC predictions 



The standard physical scale approach is to set M = = Q and adopt MS subtraction. 
Setting M = /i implies that a = d, and hence the moments have the form, 

M{n;Q^) = Ani-^) {l + Ria + R2a^ + . . .) . (24) 
\ 1 + ca / 

The coefficients Ri can be determined by expanding Eq. ([8]) in powers of a, 

Ri = + ^ (25) 

dj cdi ndi dc2 d2 , . 

= ^^+26^-^ + ^-^+26' 

and the coupling in this expression is the three- loop MS coupling with fi = Q. 

The CORGI idea (see Ref. [6j for a detailed discussion) is that all RG-predictable 
information about higher perturbative coefficients, available at a given fixed-order of calcu- 
lation should be resummed to all-orders. Given an NLO calculation for instance one knows 
Xo(Q) but not X2, or higher FRS-invariants. One should therefore resum to all-orders 
all the terms not involving these unknown invariants. As discussed in Section 2 these terms 
are multinomials in ri,fi,C2, ■■■ ,Ci,di,d2, ... ,di, .. .. 

Crucially this all-orders sum must be FRS-invariant, as separately must be the subset 
of terms involving X2, X-^ .... One may exploit this invariance and choose to use the FRS 
where all the FRS parameters are zero, ri = fi = C2 = . . . = Ci = . . . = di = d2 = ■ ■ ■ = di = 
. . . = 0. 

Setting ri = fi = means that ^ = M, setting di = then implies that (from 
Eq. (fT9l) ) tm = 61n((5/A_A/(). Also, with q = dj = the integral X(a) of Eq. ([9]) vanishes, and 
one finally obtains the CORGI form of the moments, 

M{n;Q^) = A„ (l + Xsag + Xgag + . . .) . (27) 

V i + cao / 



Here the CORGI coupling is the coupling in a 't Hooft scheme [T^, in which q = (i > 1). 
This can be written in terms of the Lambert W function defined implicitly by W{z)e^^^^ = z 

mm, 



with. 
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n 


^CORGI 


X2{n) 




n 




X2{n) 


1 


0.4688 


-1 




11 


0.2749 


-3.776 


2 


0.7156 


-3.01 




12 


0.2659 


-3.749 


3 


0.5074 


-3.28 




13 


0.2580 


-3.716 


4 


0.42966 


-3.485 




14 


0.2510 


-3.677 


5 


0.3838 


-3.627 




15 


0.2447 


-3.633 


6 


0.3530 


-3.713 




16 


0.2390 


-3.586 


7 


0.3300 


-3.766 




17 


0.2338 


-3.536 


8 


0.312 


-3.792 




18 


0.2291 


-3.483 


9 


0.2974 


-3.8 




19 


0.2248 


-3.428 


10 


0.2853 


-3.793 




20 


0.2207 


-3.372 



Table 1: The numerical values of a;coRGi and the NNLO CORGI invariants X2{n) for the 
n = 1 — 20 moments of F^. 

W-i refers to the branch of the Lambert W function required for asymptotic freedom, the 
nomenclature being that of Ref. [20]. Aj^ is the invariant scale connected with the Xq{Q) 
FRS-invariant, defined in Eq. (jl9p . Since it is an FRS-invariant it can be evaluated in any 
FRS. Choosing the MS scheme with M = = Q one finds 

Am = Am.(|) '''e-p{| + 7}, (30) 

with ri and di calculated in MS with M = fj, = Q. The factor of (^) ^^'^ converts to 
the standard convention for integrating the beta-function equation and defining Ams (see 
Ref. [21] for further details). The second factor on the r.h.s. of Eq. (|27p resums to all-orders 
the RG-predictable terms not involving X2,X3, .... The Og term sums to all-orders the RG- 
predictable terms involving X2, but not X3,X4, . . ., etc. 

The CORGI result corresponds to an MS scale choice M = fi = xQ, with 

X = Xcokg: - «-p{-^-7}- (31) 

To illustrate how the CORGI scale differs from the PS choice {x = 1) we plot in table 1 the 
XcoRGi) for the first 20 moments n = 1, 2, . . . , 20. We also tabulate the corresponding X2{n) 
NNLO CORGI invariants obtained from Eq. (j23p . The anomalous dimension coefficients up 
to NNLO are taken from Refs. [3 [9], and the coefficient function from Ref. [7]. We assume 
Nf = 5 active quark flavours. We see from table [T] that as n increases the CORGI scale 
decreases, becoming significantly less than x = 1 (PS). The X2{n) invariants are seen to be 
moderate in size. 

Finally, we discuss the third variant of perturbative QCD which we shall consider. By 
setting M = fi and rearranging, we can recast the perturbation series for M{n] Q^) of Eq. ([8|) 
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in the form 



\ d/h 

M{n;Q') = An[cR{a)) , (32) 



where 

R{a) = a + Ria"^ + R2a^ , (33) 

is an effective charge |10] . and the coefficients Ri have the form, 

~ _ bcRi 2 

Ki — — c 

a 

bcr-i (Itc n , 



bcR2 b^Rjc bRjc bRi 



c 



„3 



d ^ 2cP 2d d 
bcr2 bdiTic d2C C2C 

^^r^c ^r^c d^c bric^ SdiC^ 3 



Here i?i and R2 are the coefficients defined in Eqs. (j25p and (j26p . Rather than integrating the 
effective charge beta-function we shall instead apply CORGI to the effective charge, avoiding 
the need to numerically solve a transcendental equation which would make the fitting to data 
considerably more complicated. At NLO the CORGI and EC results agree exactly. We have 
the CORGI result 

M(n;Q2) = Ay'"' {ao + X2al + X^a^ + . . (36) 

In this case, the Xi coefficients are the CORGI invariants corresponding to single scale RS- 
dependence [6]. They have the form, 

X2 = R2 - Rj - cRi + C2, (37) 
X3 = R3-3RiR2 + 2Rl + ^-RiC2 + -C3. (38) 

The CORGI coupling oq is that of Eq. but with the scale Am now defined by, 

= (^) '^'expf^^Ms. (39) 



V b J " \^ b ^ 

We shall refer to this variant of perturbation theory as 'EC for simplicity, even though as 
noted above it is really CORGI applied to a single-scale effective charge. 

We note that we can streamline the calculation of the FRS-invariants Xi by using 
the single-scale effective charge. If we set M = /i in Eq. ([8]), then the moments reduce to 
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a single-scale problem [6]. We can then rearrange the resultant expression in terms of an 
effective charge, R{a), 

MM') = aA-^^^] . (40) 

R{a) has the form, 

R{a) = a + Ria' + R2a^ + R^a'^ + .... (41) 

The coefficients Ri can be determined by expanding Eqs. (f24l) and (fiOl) in powers of a and 
then equating coefficients. They are found to be, 

Ri = ^Ri, (42) 
a 



(43) 

where Ri and are given by Eqs. (j25p and (j26|) . If we then CORGI-ize this effective charge, 
we have a new set of FRS invariants [S] , 

M 

Xo = blnj--Ri, (44) 
^2 = R2 - Rj - cRi + C2, (45) 



and the moments become. 



d/b 



Min;Q') = { ^— ^1 . (46) 

1 + c ao + X2af, 



Expanding this into a form which we can compare with Eq. (j27p . gives, 

cat) \'^^^ f ^ , d ^ 2 



M{n;Q') = An{jf^] il + ^X2at + ...]. (47) 



Isolating the ©(cg) in the RHS bracket of the above equation, and then using Eqs. (HSI) . (H2 
(03]), (USD and (IMD, gives, 



^X2 = ^[R2-RI-cRi + C2) (48) 

(1 ^ \ 9 di dc2 (i? cdi d2 , . . 

So we see that the coefficient of the 0(0^) term in Eq. (j47p is the FRS invariant X2, of 
Eq. (|23p with fi = M (ri = n). Isolating the term will yield X3, and so on for higher Xi. 
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The coupling in Eq. (j47p is the 't Hooft couphng of Eq. (|28|) . but with the scale 
parameter determined by Eq. (j44p . Evaluating Eq. (j44p in MS with M = = Q and using 
the standard definition of the single-scale RS invariant Xq [101 [14] gives, 

Xo(Q) = 61n-^ (50) 

= 61n2-i?i. (51) 
A 

Comparison with Eq. (fT9]) then reveals that Xo((5) = (h/ d)XQ{Q). Using the same procedure 
we can obtain expressions for X3 and higher CORGI invariants. 



Non-perturbative effects 

The three variants of NNLO perturbative QCD, PS, CORGI, and EC, can all be computed 
given MS anomalous dimension coefficients up to NNLO [HI [9], and the coefficient function 
[?]• However, these perturbative predictions will be subject to non-perturbative corrections 
in the form of O (l/Q^) terms. The two principal sources of these terms are: higher twist 
terms and effects due to the mass of the target hadron. 

The perturbative form of the moments is derived under the assumption that the 
mass of the target hadron is zero (in the limit 00). At intermediate and low 

this assumption will begin to break down and the moments will be subject to potentially 
significant power corrections, of order O {rnf^/Q"^^, where mj\i is the mass of the nucleon. 
These are known as target mass corrections (TMCs) and when included, the F3 moments 
have the form [221123] . 

A^™-(n;Q2) = A^(n; Q^) + !^^J^^A^(n + 2; Q^) + O . (52) 

The moments will also be subject to corrections from sub-leading twist contributions 
to the OPE. These effects are poorly understood and hence we only estimate them; this is 
done by means of an unknown parameter, ^ht- The estimate has the form [3], 

>t«^(n;Q2) = ri(^A^^^^M{n;Q'), (53) 

and the value of A^'^ is obtained by fitting to data. Due to the poorly understood nature 
of these effects, we do not include the above term in the full analysis. Rather, we perform 
the analysis with and without this term included, and take the difference in the results as an 
estimate of the error associated with our ignorance of the true nature of these effects. 

The bottom quark mass threshold is within the range of spanned by the available 
data for F3. It is therefore necessary to evolve the expressions for the moments over this 
threshold, and in order to do this we use the formalism of Ref. [2S]. We use massless QCD 
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with 4 quarks for < and massless QCD with 5 quarks for Q"^ > m^. Here m;, is the 
pole mass of the 6-quark with mi, = 4.85 ± 0.15 MeV j24j . From the decouphng theorem, one 
finds the foHowing relation between the coupling above and below a quark threshold (denoted 
by af^i{Q'^) and af{Q'^) respectively) [25] 



af{ml) = af+i{ml) + — {af+i{ml)) . (54) 

In practice, this matching is implemented by adopting different values of the scale parameter 
in different Nj regions. This is governed by the following equations 

= ^N,[^^] x^^P( 26^,+! j' (55) 

where 5i 




4(c^/+i - c^/ ) In - 4c^/+i In (56) 
PVT^V^- ^ ^lnL.+ (c^/+^)^-(c^/)^ 

. cf^-V.^). (-) 

Here, Aat^. is the scale parameter in the region where Nj quarks are active, rriAr^. is the pole 
mass of the / quark, b^f , c^f and c^-^ are simply b, c and C2 evaluated for A'^^ quark flavours 
and we have defined = In ^m^^^^/A^^^ . Furthermore, we also demand continuity of the 
moments at the threshold i.e. 

M{n;ml)\^^,^^ = M{n;ml)\j^^^^ . (58) 

As a consequence of this, the parameters An also have different values in the Nf = 4 and 
Nf = 5 regions and their values are related by, 



/l(5) 



An 



M{n;ml) 



Nf 




M{n;ml) 



An 



^^='. (59) 



TV, 



4 The method of Bernstein averages 

When comparing theoretical predictions for moments of structure functions with experimental 
data, we are faced with the long-standing issue of missing data regions at high and low x for 
low and high respectively. This is demonstrated in Fig. [1] in which we plot the CCFR 
data [l] for 12 different values of Q^. We can see that at the lower range of we are limited 
to low-a; data, and that at high we are limited to the high x range. 

In order to reliably evaluate a moment at a particular Q^, we require data for the 
whole range of x. This being unavailable, we are forced to make some guess about how the 
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Figure 1: Data for xF^ plotted against x for the 12 different bins of the CCFR data. 
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structure function behaves in the missing data region. That is to say, we have to choose some 
method of modelhng (extrapolating and interpolating) the data to cover the full range of x. 
However we wish to make the evaluation of the experimental moments as free from QCD 
input as possible, thus making the comparison between theory and experiment as direct as 
possible. To this end, we shall adopt the approach involving Bernstein averages [21 H]; objects 
which, though related to the moments, have negligible dependence on the modelling method 
adopted (and hence on the behaviour of the structure function in the missing data regions). 

We define the Bernstein polynomials as follows, 

Pnk{x^) = 2:p7— tx''(1-^Y-', n,fcGL (60) 
r (fc + r (n — A; + 1) 

These functions are constructed such that they are zero at the endpoints x = and x = 1, 
and they are also normalized such that pnk{x)dx = 1. Furthermore, if we constrain n 
and k such that n > k > 0, then Pnk{x) are peaked sharply in some region between the two 
endpoints. 

The Bernstein polynomials can be treated as a distribution, with a mean, 

Xnk = [ XPnk{x)dx (61) 

(62) 







r(A: + l)r(n + I) 



r(A: + i)r(n + 2)' 
and variance, 

AXnk = / {X - Xnkf Pnk{x) dx (63) 
^0 



k + ^ /r(A: + i)r(n + |; 



n+l \T{k + ^)T{n + 2) 



(64) 



The Bernstein averages of -F3 are then defined by. 



FnkiQ^) = [ Pnk(.x^)F3{x,Q^)dx. (65) 
^0 

Thus F„fc(Q2) is the average of the structure function weighted such that the region around 
Xnk is emphasized. By picking the values of n and k wisely, we can construct a set of averages 
which enhance the region for which we have data for F3 and de-emphasize the regions where 
there are gaps. Therefore, in the resultant averages, the dependence on the missing data 
regions will be heavily suppressed. 

Defining this more carefully, for a given value of Q^, we only consider averages for 
which the range. 



Xnk 



\l AXnk < V'^Xnk, (66) 
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Figure 2: Constructing the Bernstein average, F^2{Q'^ = 50.1 GeV^). The hght grey region 
represents the interval in Eq. (I66p and the dark grey areas represent the missing data regions. 
The small size of the dark grey region in the right hand plot demonstrates that this average 
will have negligible dependence on the missing data regions. Note that the right hand plot 
actually shows the integrand of the Bernstein average. The average itself will be this function 
integrated over [0, 1]. 



lies entirely within the region for which we have data. The only exception to this is that if 
the highest-x data point lies within this range, then we do accept this average, but only if 
the data suggests that xFj, vanishes rapidly beyond this point. 

The construction of an acceptable average, and the resultant suppression of the missing 
data region is demonstrated in Fig.O We see that the shaded (dark grey) missing data regions 
almost disappear in the right hand plot. 

By expanding the integrand of Eq. (j65p in powers of x, and using Eq. ([T]), we can 
relate the averages directly to the moments. 



and so theoretical predictions for the averages can be obtained by substitution of Eqs. (p< 
()24p and (j36p into the above expression. The Bernstein average is seen to be a linear combi- 
nation of odd moments. Due to the unavailability of results for d2 for even n, previous NNLO 
analyses of this kind have been limited to the inclusion of only odd moments. However, 
now that the NNLO calculation of the NS anomalous dimension is complete [SJ [S] , we are no 
longer constrained in such a way. In light of this, we define a new set of modified Bernstein 
polynomials, 

which include only odd powers of x and hence whose averages are related to even moments. 



These modified Bernstein polynomials are simply the original polynomials of Eq. (j60p . mul- 
tiplied be X, and then 're-normalized' such that they still satisfy Jq pnk{x)dx = 1. We can 
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calculate the mean and variance of Pnki^)^ 



_ r(fc + |)r(n + 2) 
" r(fc + i)r(n + |)' 



Ax 



k + 1 
n + 2 



r(A; + |)r(n + 2) 



(70) 



r(A;+l)r(n + § 

and in analogy with Eq. (j65p . we define the modified Bernstein averages, 

Fnfe(Q') = [' Pnk{x^)F3{x,Q^)dx. (71) 

Jo 

Again, we only accept experimental modified averages for which the range, 

Xnk - \/AXnk < X < Xnk + \/ AXnk, (72) 

lies within the region for which we have data. We obtain theoretical predictions for the 
modified averages using the equation, 

which is seen to be a linear combination of even moments. 

In order to calculate averages from data for F3, we need an expression for X-F3 covering 
the entire range of x, for each value of Q^. As mentioned above, the values of the moments 
calculated in this way will depend on how we model the structure functions in the missing 
data regions but for the averages, this dependence is suppressed. However, we would like to 
test this assertion, and so we use four different methods of modelling and perform our 
analysis separately for each method. Significant differences between the results would signify 
a failure of the Bernstein average method, and in instances where this is the case, we reject 
that particular average at that particular Q^. Moderate deviation however, is acceptable, 
provided that we use the magnitude of the deviation as an estimate of the error associated 
with the missing data region. This error is then included as a 'modelling error' in the final 
result. In this way, we can almost completely remove any dependence on missing data regions, 
and also quantify the error associated with any residual dependence. 

The four extrapolation methods we use are described below: 

I In the first method, we fit the function, 

xF^ix) = Ax'^{l-xf, (74) 

to the data for each fixed value of Q^. The parameters B and C are obtained by 
performing fitting of Eq. (I74p to data for F^. They are Q^-dependent quantities, 
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Figure 3: Diagram depicting the x-ranges covered by the CCFR data, at different Q^. 

and errors on their values are obtained by performing the fitting with the data for F3 
shifted to the two extremes of the error bars. 

A justification for the particular form of fitting function in Eq. ()74p can be found in 
Ref. [30]. However, the simple fact that this function fits the data well is justification 
enough, since the Bernstein averages are independent of the extrapolation method. 

II The second method we use is linear interpolation between successive data points. We 
also extrapolate beyond the data range, to the endpoints xF3{x)\x=o = xF3{x)\x=i = 0, 
in order to be consistent with method I. 

III The third method consists of using the fitting function of Eq. ()74p . but setting xF^^x) = 
everywhere outside the region for which we have data. 

IV In analogy with III, in this method we use the linear interpolation of method II but 
setting xF3{x) = everywhere outside the data region. 

The deviation between the results obtained from the above methods (in particular, the 
diff'erence between the first two and the last two) will be a good measure of the effectiveness 
of the Bernstein average method. 

Data for xF^ in neutrino-nucleon scattering is available from the CCFR collaboration 
[T] . The data was obtained from the scattering of neutrinos off iron nuclei and the measure- 
ments span the ranges 1.26 < < 199.5 GeV^ and 0.015 < x < 0.75. The x-ranges covered 
at each are depicted in Fig. [3l 

In Fig. H] we show each of the four modelling methods applied to F^ measured at 
= 79.4 GeV^. Also shown on these figures (in grey) are the fits which are used to 
determine the errors on the modelling, which propagate through to errors on the averages. 
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Figure 4: The four methods used for fitting the structure functions. Here we show the 
measured values of xF^ at = 79.4 GeV^. The errors are determined by re-performing the 
fitting for the data shifted to the extremes of the error bars, and this is denoted by grey lines. 



From the CCFR (and using the methods I - IV outlined above) we can obtain expres- 
sions describing the behaviour of the structure function over the full range of x, for each value 
of Q^. It is then possible to extract experimental values of the averages, using the methods 
outlined below: 



In the case of I, obtaining the averages is particularly simple. Substituting Eq. (I74p 
into Eqs. ([65]) and ([7T|) gives, 

= ^f|^E«^;i^B(2(..0 + fi.C + l), (75) 



for the Bernstein averages and. 



-(exp) ^ 2r(n + 2) ^' (-1) ' 

' 1) ^ lUn-k-iy. 

' 1=0 ^ ' 



for the modified Bernstein averages. Here, B(x, y) = Y {x)Y {%)) jV [x + y) is the Beta function. 
Once values for A, B, and C, have been obtained, substitution into the above expressions 
leads directly to the averages. 



18 



In the case of II, each of the averages is spht into j + 1 sections (where j is the number 
of data points) and each section is an integral over a polynomial of order 2n + 1. It is then 
reasonably simple to evaluate the averages by computing this set of integrals. This approach 
also applies to method IV, but in this case there are only j — 1 integrals. 

For method III we simply integrate the fitting function, multiplied by the Bernstein 
polynomials, with the integration limits being the values of x at the first and last data point. 

Having outlined the method for obtaining the experimental averages we now turn our 
attention to which averages are acceptable at which energies. The highest moment we use 
is the 18th moment and the lowest the 1st. Inclusion of higher moments than this leads 
to no significant increase in the number of acceptable Bernstein averages. The upper limit 
of n = 18 implies that the highest Bernstein averages included are Fgi^ and -Fgfc ^-rid that 
the lowest used are Fio and Fiq. We exclude the averages for which n = k as they simply 
correspond to individual moments themselves. This leaves us with a total of 72 (36 + 36) 
potential averages at our disposal for each value of Q^. This number will be reduced when 
we come to exclude averages on the basis of the acceptance criteria. After applying the 
acceptance criteria, we are left with 132 data points for the standard Bernstein averages and 
141 for the modified Bernstein averages. Exactly which averages we use at a particular Q^, 
can be determined by inspecting the plots in the results section. 

In Fig. [5] we plot the dominant regions of the Bernstein polynomials (given by Eq. (I66p ) 
for each of the used averages. This is superimposed onto the data-range diagram of Fig. [3l 
Figure [6] shows equivalent plots for the modified Bernstein averages. These plots can be used 
to identify which averages are acceptable for a particular value of Q^. 

5 Fitting procedure 

We use minimization to optimize the fits of the theoretical predictions to the data. The 
highest moment included in the experimental averages is the 18th, and so when TMCs are 
included, we will require predictions for the first 20 moments. Therefore, the set of fitting 
parameters comprises of {^i . . . ^20} plus the QCD scale parameter A^g. When we include 
higher twist corrections this set is expanded to include Aht- To check consistency between the 
odd and even moments we perform the analysis for each of these sets of moments separately 
and then finally together, and compare the results. 

Although we stated previously that 132 standard Bernstein averages are available 
to us, in the CORGI case this is reduced to 130 for the following reason: When fitting 
predictions to the data, we scan values of A^ between and 590 MeV for a minimum 
in x^- Unfortunately, for n = 17 and 19, the values of A_a^/Ams (see table [1]) are such that 
= 7.9 GeV^ is below the Landau pole in Eq. ([281) . Consequently we cannot obtain CORGI 
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Figure 5: The black and light grey bars (I — \) show x ranges covered by the CCFR data at 
different energies. Superimposed onto these, in various colours, are the peaked regions of the 
individual Bernstein polynomials, defined by the interval in Eq. (j66p . 

predictions for Bernstein averages which include these moments. Prom Eqs. ()52p and (j67p we 
can determine that this excludes FsoiQ"^ = 7.9 GeV^) and FjoiQ"^ = 7.9 GeV^) from the fit. 
In the case of the modified Bernstein averages, Fgo 

(Q2 = 7.9 GeV^) and ^70(^2 = 7.9 QeV^) 
are already excluded due to their failure to meet the acceptance criteria. 

The CCFR data includes statistical errors and 18 different sources of systematic error. 
These errors cannot be added in quadrature, and so we perform the analysis for each of these 
19 sources of error separately and then add the variation in the results in quadrature to obtain 
the final total error. We also include, as additional sources of error the deviation in results 
associated with using the four different modelling methods (this forms the 'modelling error' 
in our final result), and the deviation in the results obtained from performing the analysis 
with and without HT corrections included (forming the 'HT error'). 

Correlation of errors 

When fitting theoretical predictions to experimental data using minimization, care must 
be taken in order to take into account fully the correlation between data points. 
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Figure 6: The black and light grey bars (I — \) show x ranges covered by the CCFR data at 
different energies. Superimposed onto these, in various colours, are the peaked regions of the 
individual modified Bernstein polynomials, defined by the interval in Eq. ()72|) . 

To construct from a set of uncorrelated data points {Z""'' } (i = 1, . . . N), with 
errors and corresponding theoretical predictions {fi^""'}, we have, 

^ / f exp. ^thoo. \ 2 

= E( ^ - ' ) ■ (77) 



i=l 



The raw data for xF^ are uncorrelated. However, we are not comparing predictions for the 
structure functions themselves with data directly; rather we are doing so indirectly via the 
Bernstein averages. For a given value of Q^, the full set of Bernstein averages (and modified 
Bernstein averages) we obtain will be correlated, due to their being derived from the same 
set of {xF^) data points (see Fig. [2]). 

In the case where {/™'' } are correlated the function becomes [26], 

N N 



i=i j=i 

^^exp. ^thco.'^^ ■^■—1 ^^cxp. ^thco.^ ^78^ 
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In the second line of the above equation we have constructed vectors from the data points and 
their predictions. V is known as the covariance matrix. It encodes the correlation between 
each of the data points; its elements are obtained as follows, 

Vij = cov {fi, fj) 

= {m-{n){fj)- (79) 

If the fi are themselves functions of M variables (representing the 'raw' data), then we 
have, 

cov (A,/,) = ||(||)(^)„v(..,x,). (80) 

If the data for Xi are uncorrelated, this reduces to. 




To obtain the covariance of the Bernstein averages we simply substitute -Fnfc(Q^) for fi and 
F:i[x) for Xi in Eq. (jSip . We assume no correlation between different values of and so 
Eq. ([75]) decouples into 7 different matrix equations (one for each of the values of between 
7.9 and 125.9 GeV^). By using the trapezium rule to approximate the Bernstein average 
integrals, we obtain the following expression for the covariance matrices: 

1 

VlrniQl) = XI 4 (Pnkixj))i {Pnk{Xj))^ [Xj+l - Xj^if CF^^^^^y (82) 

Here, the index j runs over the number of data points we have for xF^ at a given Q?. The 
J = and j = Ni + 1 terms are simply the x = and x = 1 endpoints (see parts I and II 
of Fig. S]) . From this equation we can calculate elements of the covariance matrices for each 
value of Q^. All that remains is for us to invert them. However, upon attempting to do so, 
we find that these matrices are ill-conditioned, with some of their eigenvalues being close to 
zero. Hence their inverses are intractable. As a result of this, it is impossible to perform a 
reliable analysis of the averages with their correlation taken into account. 

We believe that this is principally due to the fact that the correlation between averages 
is significant in some cases, and this in turn is an artefact of the fact that the selection criteria 
systematically select Bernstein polynomials which are peaked in the same region and hence 
are of fairly similar shape. This situation arises because the intent behind the inclusion of 
more averages in the analysis is not to increase the amount of 'data', rather it is to further 
ensure that the missing data regions are suppressed. 

In light of this, we settle for the method adopted in Refs. (Sj S], in which the naive 
function of Eq. ()77p is used, but the error bars on the averages are modified in order 



22 



to account for the 'over-counting of degrees of freedom'. For example, for the standard 
averages, at each value of we have theoretical information on 9 moments, but the number 
of experimental Bernstein averages we use is often more than this; e.g. for = 20 GeV^ 
we use 27 standard Bernstein averages. To remedy this, we adopt the following approach. 
For each value of we count the number of averages above 9 as duplicate information. 
The number of duplicates we have altogether is 73 and so for the values of for which 
we have more Bernstein averages than moments, we rescale the error on these averages by 
y^l30/ (130 — 73) = 1.510. Correspondingly, for the modified Bernstein averages we rescale 
the errors by a factor of y^141/(141 — 87) = 1.616. This rescaling has the effect of suppressing 
the contribution of the duplicate data points to relative to those values of for which 
we have fewer Bernstein averages than moments. 



Positivity constraints 

The fact that X-F3 is a positive definite function, and that the moments are simply integrals 
over these functions multiplied by a single power of x, means that we can impose certain 
positivity constraints on the parameters as follows. 

We construct the following matrices from the moments, 

I Mx M2 ■■■ M9 \ 

M = M3 ■•• , (83) 

V Mg Mn I 

and 

/ A>Ii A7W2 ••• AA^9 \ 

A7W = ■■• , (84) 

y ATWg A7Wi7 j 

where Mn = M{n; Q"^) and AMn = Mn - Mn+i- 

In order for M{n;Q'^) to be moments of positive definite functions (as the structure 
functions must be), the determinants of the above matrices, and of all their minors, must 
be positive, for all values of [4]. Evaluating these determinants at fixed will translate 
to conditions on the parameters An- We do not implement these constraints as part of the 
fitting procedure. Rather, we perform checks on the values of the fitting parameters resulting 
from the minimization to ensure that they obey the above constraints. 
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However, we do impose positivity constraints on the moments themselves. As a result 
of the determinantal constraints described above, and from the general form of the moments 
given in Eq. ([T]), we can infer that the following inequalities must be satisfied, 

Min;Q'^) > 0, (85) 
7W(n;Q2) ^ Min + 1;Q^), (86) 

for fixed Q^. Furthermore, we can implement these constraints by defining our fitting param- 
eters An in terms of a new set of parameters and then minimizing with respect to these 
new parameters. 

We begin by picking some value of Qq at which to implement the conditions. We then 
take the last moment used in the analysis (n = 20) and rewrite the constraint in Eq. (j85p as, 

M{20;Ql) = (iso)', (87) 

where A20 is a real number. The constraints in Eq. (186p can also be rewritten as, 

M{n;Ql) = M{n + l;Ql) + (^A^y , (88) 

for 1 < n < 20, where An are all real numbers. The LHSs of Eqs. (1870 and (I88p are simply 
a fitting parameter times a number. For example, in the case of n = 2 and Qq = 12.6 GeV^ 
we have, 

A^(2; 12.6 GeV^) = 0.3932^2. (89) 

Prom this, and equivalent expressions for the rest of An, we can obtain an expression for each 
An in terms of the parameters Ai - ^20- This means that we can replace the parameters Ai - 
A20 with Ai - A20 in the function. By doing this and then minimizing with respect to the 
An parameters, we can find a minimum in for which the constraints in Eqs. (j85p and ()86p 
are automatically satisfied. In effect, the reparameterization embedded in Eqs. (f87l) and (f88]) 
restricts the parameter space to exclude solutions for which the constraints are not satisfied. 

To implement this reparameterization we must choose a value of at which to 
impose the constraints, whereas in reality they must be satisfied for all Q^. Because of this, 
we perform the analysis for several different values of Qq and check that the results remain 
stable. 



6 Results of fitting to the data 

We focus principally on the results from the CORGI analysis in which both odd and even 
moments are included and in which we include target mass corrections. This analysis results 
in a prediction for the QCD scale parameter of, 

Ag = 219.11^22.1 MeV. (90) 
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The errors on this value can be broken down into four different sources, 

Ag = 219.11 ™ (,tat.) i|f7(sys.) t\tf, (mod.) ± 8.97 (HT) MeV. (91) 

We have used method I to obtain experimental values of the averages. The deviation between 
the results obtained using methods I and IV is used to evaluate the modelling error since 
these are the two methods which exhibit the largest deviation. 

This result for A^g corresponds to a value of the strong coupling constant (evaluated 
at the mass of the Z particle) of, 

as{Mz) = O.nSdtSl (92) 

These values are in excellent agreement with the current global averages of A^ = 207.2 GeV 
and as = 0.118 it .002 [23|. There is also good agreement with the result obtained from 
fits using the Jacobi polynomial method [13] which yield as{Mz) = O.IIQIqqq^. This result 
is based on fits using odd moments only and includes a contribution to the error from scale 
dependence. Whilst close to the global average, the value of as{Mz) we obtain is significantly 
larger than the values found in a recent analysis of the F2(x,Q^) structure function [27], or 
from fits of parton distibution functions where DIS and Drell-Yan data are combined [28], it 
is also larger than the value which minimizes the global parton distribution function 

fits such as Ref. [29] . 

The x^/d.o.f. for our CORGI result is as follows, 

^ 20.37 
d.o.f. 271 -(20 + 1) 

= 0.0815. (93) 

Here '271' refers to the number of experimental Bernstein average points used in the fits. 
Although this value is an order of magnitude larger than the x^/d.o.f. obtained in Ref. [1], 
it is still significantly smaller than one would expect, suggesting the errors on the Bernstein 
averages have been over estimated. However, as discussed previously, the function we are 
using does not take into account the correlation between data points. Indeed, if correlation 
was taken into account, one might expect that a more reasonable value of would be 
obtained. 

Furthermore, in the function we eventually used, the errors on the Bernstein av- 
erages were rescaled in order to take into account the 'over-counting of degrees of freedom'. 
As a result, the number of Bernstein averages is not representative of the true number of 
degrees of freedom in this particular x^ function. Indeed, the Bernstein averages in the plots 
in Fig. [7] can be constructed from just 58 different moments at different values of (via 
Eq. (j67p ). Similarly, the modified Bernstein averages in Fig. [8] can be built from 53 different 
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,2 / J „ f 

X /d.o.t. 


All moments 


219.1+1:;^ 


"•11^9-0.0019 


20.37/(271 - (20 + 1)) 


Odd moments 


210.5+|;5 


n 1 1QO+0.0029 
"•11^^-0.0030 


10.94/(130 - (10 + 1)) 


Even moments 




"•119« -0.0048 


9.24/(141 - (10 + 1)) 


All moments: > only 


232,4+34.9 

-33.2 


1200+°-°°2^ 

'^•-^^'^"-0.0027 


15.59/(228 - (20 + 1)) 



Table 2: In this table we present the result of the analysis performed using the CORGI 
approach to perturbation theory with target mass corrections included. We compare the 
results obtained when we include all moments (up to n = 20) with those obtained when we 
restrict the analysis to even or odd moments only. We also show the results from performing 
the analysis with only data points for which > ni^ {Nf = 5) included. 



moments. Hence, 

20.37 



0.226, (94) 



d.o.f. 111-21 

is more representative of the true value of x^/d.o.f. in this approach. This is a more acceptable 
value, however we stress that the true minimum in can only be determined by taking 
correlation fully into account. 

In Fig. [7] we plot the CORGI predictions for the Bernstein averages (with TMCs 
included) fitted to the experimental values. Figure [8] shows equivalent plots for the modified 
averages. 

In table[2]we present the full set of results from the CORGI analysis. This table shows 
the results obtained when we include both odd and even moments (standard and modified 
Bernstein averages) together and also when we restrict the analysis to odd and even or odd 
moments only. These results allow us to check consistency between the odd, even and 'All' 
analyses. 

We can also use the results of the 'odd moments' analysis to check consistency with 
previous analyses. In the PS analysis of Ref. (in which only the n = 1,3,5,7,9,11 and 
13 moments were included) a value of A^ = 255 ± 72MeV was found, corresponding to 
A^ = 178t55MeV, with a value of x^/d.o.f. = 0.007. Using the same set of moments, the 
CORGI analysis of Ref. [5] found a value of A^^ = 228+^^MeV, although it must be noted 
that this analysis used incorrect values of the coefficients X2 , and therefore the result must 
be regarded as unreliable. Both of those results are indeed consistent with the 'odd moments' 
analysis performed here. 

We also perform an analysis in which we restrict the CCFR data to > only, as a 
check on our method of evolving through the h quark threshold. These results are also included 
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Figure 7: CORGI fits for the Bernstein averages, with TMCs included. 



27 




Figure 8: CORGI fits for the modified Bernstein averages, with TMCs included. 

in table [21 In table 3 we give the fitted CORGI values for the An non-perturbative coefficients 
for n = 1 — 20, together with the values of the corresponding moments at = 8.75 GeV^ 
and 12.6 GeV^. 

In table H] we compare the CORGI results with those obtained using the PS and 
EC approaches. We also present results obtained from performing these analyses with and 
without target mass corrections. The fact that the number of d.o.f. for the CORGI fits is 
271 (272), as opposed to 273 for PS and EC, refiects the fact that for the smallest energy 
bin Q2 = 7.9 GeV^, the appearing in the CORGI coupling exceeds for the highest 
n = 19, 20 moments, and hence one is below the Landau pole in the CORGI coupling of 
Eq. (p8]) . Correspondingly Xcorgi is significantly less than unity (see table 1). We simply 
omit the two affected Bernstein average points from the CORGI fit. 
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n 


4(4) 


>I(n;8.75 GeV^) 


M{n; 12.6 GeV^) 


1 


2.346 


2.494 


2.525 


2 


0.8814 


0.3557 


0.3466 


3 


0.4133 


0.1002 


9.545x10-^ 


4 


0.2217 


3.835x10-^ 


3.584x10"^ 


5 


0.1292 


1.744x10"^ 


1.603x10-^ 


6 


8.134x10"^ 


9.048x10"^ 


8.191x10-^ 


7 


5.241x10"^ 


4.988x10"^ 


4.452x10-^ 


8 


3.639x10"^ 


3.044x10"^ 


2.681x10-^ 


9 


2.434x10"^ 


1.826x10"^ 


1.588x10-^ 


10 


1.822x10"^ 


1.246x10"^ 


1.07x10-^ 


11 


1.202x10"^ 


7.588x10^4 


6.438x10-4 


12 


9.64x10"^ 


5.677x10-4 


4.76x10-4 


13 


5.935x10"^ 


3.289x10^4 


2.726x10-4 


14 


5.119x10^^ 


2.69x10-4 


2.204x10-4 


15 


2.702x10"^ 


1.355x10-4 


1.098x10-4 


16 


2.535x10"^ 


1.22xl0"4 


9.768x10-^ 


17 


9.362x10""^ 


4.343x10-'' 


3.44x10-*^ 


18 


9.739x10-'^ 


4.376x10-^ 


3.426x10"^ 


19 


9.807x10-^" 


4.284x10"" 


3.317x10"" 


20 


7.69x10"^" 


3.277x10-" 


2.509x10-" 



Table 3: Fitted values of An (in the Nf = 4 region) together with the moments evaluated at 
Q2 = 8.75 and 12.6 GeV^, for the CORGI approach. 





AS(MeV) 


as{Mz) 


xVd.o.f. 


CORGI 


with TMC 


219.1+g; 




20.37/(271 - (20 + 1)) 


no TMC 


280.3^-:^ 


n lOQK+o.ooi? 

^•^^"^^-0.0017 


24.76/(272 - (18 + 1)) 


PS 


with TMC 


?no 4+25-8 


n 1170+0.0023 


21.73/(273 - (20 + 1)) 


no TMC 




1219+0-0020 


25.20/(273 - (18 + 1)) 


EC 


with TMC 




+0.0017 
-0.0017 


22.71/(273 - (20 + 1)) 


no TMC 


261.5;^;- 


r| -|r,r,r, +0.0014 

U--L^^^_0.0014 


26.04/(273 - (18 + 1)) 



Table 4: In this table we compare the results of the analysis performed with the three different 
approaches to perturbation theory described in section 3, CORGI, PS and EC. We also show 
the results from these analyses performed with and without target mass corrections. 
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7 Discussion and conclusions 



In this paper we have used three different approaches to perturbation theory to perform a 
phenomenological analysis of moments of using the method of Bernstein averages. The 
three approaches differ in how they deal with the FRS dependence. In the CORGI approach, 
we allow the FRS invariant quantity Xq{Q) to determine the relationship between M, fi and Q 
for each moment. In so doing, we automatically resum the subset of terms present in the full 
perturbative expansion which are RG-predictable at NNLO. In the physical scale approach 
we set M = fj, = Q and adopt the MS scheme for the subtractions in the renormalization 
and factorization procedures. In the effective charge approach, we set M = fj, and apply the 
CORGI approach to the resulting single-scale effective charge. We described how predictions 
are derived in these three approaches and corrected errors in the CORGI method which were 
present in Refs. [5l[6]. 

We described how target mass and higher twist corrections affect these theoretical 
predictions and also how we evolve expressions for the moments through the 6-quark thresh- 
old. We explained how the Bernstein averages method eliminates any potential dependence 
of the analysis on missing data regions in x and Q^, and we also described how this method 
is generalized to treat both odd and even moments. We described the fitting procedure used 
to extract the optimal values of the QCD scale parameter and how we can implement various 
constraints which ensure that the results of this fitting are consistent with the structure func- 
tions being positive definite functions. We also presented an alternative, and slightly easier 
method for deriving the FRS invariant quantities Xj. 

The results of the CORGI analysis presented in table [2] show excellent agreement with 
the current global average for the strong coupling evaluated at = |24] . and are also 
in good agreement with fits based on the Jacobi polynomial approach [13]. From this we 
conclude that CORGI perturbation theory performs well when applied to the analysis of 
moments. The analyses in which we include only odd or even moments are consistent with 
each other and with the full (all moments) analysis. Furthermore, in the analysis in which 
we include all moments, the errors are greatly reduced. This improvement in the analysis 
is made possible by the availability of the full NNLO anomalous dimension calculation and 
represents significant improvement on previous analyses. 

Excluding data points for which < leads to no significant change in the results 
and from this we conclude that the quark mass threshold method we have applied is suitable 
to the moment analysis. The error associated with the exclusion of higher-twist effects, given 
in Eq. (j9ip . is relatively small, signifying that these effects are not particularly important at 
scales Q2 > 7.6 GeV^. 

We include in the analysis positivity constraints on the moments (Eqs. (j85p and (j86p ). 
via the parameter redefinitions defined in section 4. We find that this implementation has 
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little effect on the prediction of A^g (~ 10 Mev), but does make a difference to the values of 

The CORGI predictions for the Bernstein averages (with TMCs included) are plotted 
in Figs. [7] and [8] and show excellent agreement with experimental values. This is reflected 
by the low value of x^/d.o.f. associated with this fitting, given in Eq. (I93p . However, as we 
noted, ideally the full covariance matrix should be used in constructing to account for 
correlations between the Bernstein averages used in the fits. Unfortunately, we found that 
this matrix is ill-conditioned, having some eigenvalues close to zero, and so it proved to be 
numerically intractable to invert the matrix to construct the true x^- We therefore resorted 
to the same approximate rescaling of errors employed in Refs.OHj to try to compensate for 
possible correlations. 

The results also show consistency between CORGI, PS and EC. The PS and EC 
analyses lead to values of A^g and as slightly lower than in the CORGI analysis. However, 
this variation is well within the error bars on the associated quantities. Inclusion of HT 
effects generally results in a small shift in Ams of about 10 MeV. However, when target mass 
corrections are included, we see a shift of approximately 60 MeV in the predicted value of 
Ajjs and from this we conclude that these contributions are significant in the case of Fj,. 

An obvious further study would be to apply the same fitting procedure to the recently 
released NuTev data [2]. In future work we also hope to report on similar fits to data for 
the F2 structure function [31]. This analysis is considerably more complicated due to the 
presence of an additional singlet component. 

Note Added in Proof 

At around the same time as our paper was completed a related analysis of the CCFR data 
for XF3 also employing Bernstein averages appeared [32] . 
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